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The specific heat C and the electronic and phononic thermal conductivities Ke and K,ph are 
calculated in the mixed state for magnetic fields H near Hc2 ■ The effects of supercurrent flow and 
Andreev scattering of the Abrikosov vortex lattice on the quasiparticles are taken into account. The 
resulting function C{H) is nearly linear while He{H) exhibits an upward curvature near Hc2 ■ The 
slopes decrease with impurity scattering which improves the agreement with the data on MgB2 . 
The ratio of phonon relaxation times r„/rs = giuoo, H) for phonon energy uo, which is nearly a step 
function at = 2A for the BCS state, is smeared out and tends to one for increasing H . This 
leads to a rapid reduction of Uph (H) in MgB2 for relatively small fields due to the rapid suppression 
of the smaller energy gap. 



Numerous experiments indicate that the superconducting state of MgB2 (Tc = 40K) is a conventional s-wave 
pairing state mediated by the electron-phonon interaction. Nevertheless, the theoretical explanation of physical 
quantities is demanding because two different energy gaps are formed where the larger gap Ai is associated with 
the nearly cylindrical cr-sheets of the Fermi surface and the smaller gap A2 with the three-dimensional 7r-sheets. 
Analysis of the temperature dependence of the measured thermal conductivity in the basal plane of MgB2 Q gives 
supporting evidence for two different gaps on different Fermi surface sheets. The situation is complicated because 
the thermal conductivity k is a sum of an electronic part and the phononic or lattice conductivity Kp^. In Ref. 11 
Ke and Kph were calculated with help of the BCS theory ^ where He is limited by elastic scattering (scattering by 
phonons is omitted) and Kph is limited by electron scattering. This theory for Kph is supplemented by adding to the 
phonon relaxation rate due to scattering by electrons, the relaxation rates due to scattering by point defects Q and 
other defects. We remark that the peak in Kph below Tc in the cuprate superconductors has been fitted by including 
scattering by sample boundaries, sheet-like faults, and dislocations 

Recently, the ab-plane thermal conductivity of MgB2 has been measured as a function of magnetic field H with 
orientations both parallel and perpendicular to the c-axis |6|. At low temperatures k{H) drops steeply for increasing 
_ff up to a relatively low field which we denote by {"^^ kOe), and then it rises continuously up to Hc2 (— 30 

kOe for H || c). The first drop of k is interpreted as due to the behavior of Kph{H) which is caused by a strong 
suppression of the smaller gap A2 by relatively small fields. This leads to a rapid increase in the number of those 
quasiparticles that dominate the scattering of phonons Q . The rapid suppression of the smaller gap can be explained 
in terms of a two-band model with different energy gaps where the smaller gap is induced by Cooper pair tunneling 
from the band with the larger gap 0. The observed increase of k{H) above is attributed to the increase of Ke{H) 
because Kph has already reached its normal state value |^ . Since this field dependence of k^ {H) at low temperatures 
is qualitatively similar to that of the specific heat coefhcient ^{H) [^,||, it is concluded that the growth of Ke{H) is 
due to the rapid increase of quasiparticles in the vortex cores associated with the larger energy gap. 

In this Letter we present theories for the specific heat and the electronic and phononic thermal conductivities in the 
vortex state for applied fields near the upper critical field H(,2- For simplicity we consider only one isotropic s-wave 
pairing gap either on a cylindrical Fermi surface with H || c , or on a sperical Fermi surface. Our goal is to explain the 
measured quantities ^{H) and Ke{H) near Hc2 (— 30 kOe), and Kph{H) below the effective upper critical field iJ^2 
(^ 1 kOe) for the vortex lattice associated with the smaller energy gap. Our theories are based on the normal and 
anomalous Green's functions G and F obtained from the Gorkov integral equations with kernels given by the product 
of Abrikosov vortex lattice order parameters A(ri)A*(r2) and the phase factor due to the magnetic field |]lO|. These 
Green's functions depend sensitively on svt\9 where 9 is the angle between the quasiparticle momentum p and the 
magnetic field H. For ^ 0, G and F tend to the Green's functions for a BCS superconductor while, for 9 7r/2, 
they take into account the effects of supercurrent flow and Andreev scattering due to the vortex lattice. This theory 
has been applied previously to calculate Kg (T, H) for energy gaps with line nodes like those occuring in the high-Tc 
cuprates and Sr2Ru04 Q. 

A simplified version of the theory in Ref. pO||, which has been derived from the Eilenberger equations yields 
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the following expression for the spatial average of the density of states: 
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The quantities appearing in the equation are defined as follows: z — 2[ilA + i{A/v)r]/ sm6 with 9 — Z(p,H) , 
A = {2eH)-^/'^, A = AA/z;, O = cj/A, A^ = (iJ,2 - H)/6(3aH, A^ = A2(T) [1 - (iJ/i?c2)] , and = 
exp(— 2:^)erfc(— «z) . Here v is the Fermi velocity, (3a the Abrikosov parameter which we take as 1.2, Ao(T) is the 
BCS gap, and T is the normal-state impurity scattering rate. In the Born and unitary scattering limits, T has to 
be replaced by FA or T/A, respectively, where A{n) is calculated self consistently. The specific heat C is given by 
C = NqT dx sech^(x/2)A(T.T/A) . At low temperatures (T < A), C = ls{H)T , where the specific heat coeffi- 
cient 7s (-ff) is proportional to A{Vl = 0). In Fig. 1 we have plotted A{n = 0) versus H/Hc2 for sm9 = 1 corresponding 
to a cylindrical Fermi surface and H || c, and impurity scattering rates 6 = F/Aq =0.1, 0.2, and 0.5. A is seen to 
be nearly linear near Hc2 with a slope at Hc2 that decreases for increasing S. In Fig. 1 we also show our results for 
the angular average A — J^^^ d9sm9 A{9) which corresponds to a three dimensional Fermi surface. Comparison with 
the solid curves in Fig. 1 shows that the slopes at Hc2 are decreased by the angular average. It should be pointed 
out that our results are only strictly valid in the vicinity of i?c2- However, solution of the Eilenberger equations for a 
vortex lattice shows that the spatial average of the resulting density of states is well approximated by A{^1) for fields 
down to about {l/2)Hc2 The solid curve of A vs H/Hc2 for 5 = 0.5 in Fig. 1 agrees qualitatively with the 

measured field dependence of the specific heat coefficient 7s (-ff) at low temperatures and fields near Hc2 ||,^- In the 
unitary impurity scattering limit A increases slightly for decreasing field, however at lower fields our theory does not 
apply anyway. 

We turn now to the theory of the electronic thermal conductivity Kg in the vortex state near Hc2 which has been 
developed in Ref. [0. These expressions are easily modified to apply to an isotropic s-wave pairing state and a 
field along the c-axis. Employing the expression valid at low temperatures we obtain for the ratio Kes/ ^en as a 
function of H/Hc2 the plots shown in Fig. 2 for constant sin 6* = 1 and impurity scattering rates 5 =0.1, 0.2, and 
0.5. One sees that these plots exhibit upward curvatures with decreasing slopes at Hc2 for increasing values of 6. 
In Fig. 2 we have also plotted the corresponding results for the angular averages over 9 of Kes/Ken- Comparison 
with the measured k.{H) near Hc2 , which is presumably dominated by the field dependence of Ke{H) shows that 
the measured upward curvature towards Hc2 is qualitatively best described by the upper solid curve in Fig. 2 for 
the relatively large impurity scattering rate S = 0.5. It should be noted that the renormalization of the impurity 
scattering rate in the unitary limit, F — > F/A, slightly increases the upward curvature of Kes/Ken in comparison to 
that obtained without renormalization or in the Born scattering limit. This is in somewhat better agreement with 
the data. Of course our theory cannot explain the measured field dependence of Ke{H) at lower fields which first rises 
steeply and then saturates until the upward curvature sets in at about (l/2)iJc2- Our theory also cannot describe 
the measured anomalous temperature dependence Ke oc T-^/^ for fields below approximately (1/2) Hc2 '■ i^en and thus 
Kes is proportional to T/F. 

It was argued in Ref. that the similar field dependencies of the specific heat and electronic thermal conductivity 
can be explained in terms of the relationship — C^vpilZ where I is the mean free path. In fact, our theoretical 
expression for Kg has a similar form if one sets £ = vpT^. where Te = 1/Im^o with Im^o =7 + 7^1 [0- Here 7 is the 
impurity scattering rate and 7^ is the scattering rate due to Andreev scattering of the quasiparticlcs by the vortices. 
7A is shown as a function of angle 9 for different values of and A in Ref. jlj]. These plots show that, for extended 
states (f2 > 1), 7^ increases from zero as 9 increases from to 7r/2. This means that the Andreev scattering rate 
has a maximum for quasiparticlcs moving perpendicular to the vortex axis. For increasing field (decreasing A), 
decreases which has the effect that the curve for Kes/ i^en exhibits an upward curvature in comparison to Ce which 
has a slight downward curvature for increasing field. 

Until now there exists only a phenomenological theory for the phonon heat transfer in the mixed state where the 
phonons are scattered by vortices consisting of normal state cylinders ||l^. We develop here the microscopic theory 
of Kph (limited by electron scattering) in the vortex state in close analogy to the theory for the BCS state In 
that case the sum of probabilities for absorption and emission of phonons by quasiparticlcs yields a relaxtion rate 
1/ts for a phonon of energy ujo which is proportional to the integral over quasiparticle energy E of the expression 
\EE'/ee'\[l - A^/EE'][f{E/T) - f{E'/T)]. Here E/s = E/{E^ - A^)^/^ is the density of states, E' ^ E + ujq, and 
f{E/T) is the Fermi function. In the mixed state near Hc2, the density of states E/e is replaced by A{n) as given 
in Eq. (Q). The coherence term (A/e)(A/e'), which arises from the matrix elements for quasiparticle scattering, is 
replaced by FF^ where the anomalous Green's fimction F is given in Ref. The spectral function of F yields the 
following analog of A/e: 
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The argument z was defined below Eq.(p|). The functions A{Q) and B{fl) are even and odd in Q. For 9^0 they 
tend to the BCS functions ujjiyp' — A^)"^ ^nd AjiuP' — A^)^/^, respectively. In this way we obtain for the ratio of 
phonon relaxation times in the normal and superconducting states (denoted by g in Refs. and |^): 

/>oo 

T„/r, - .9(r!o) - [1 ~ exp(-r!oA/T)] (2/l]o) / d£l /[(f! - r!o/2)A/T] /[-(f! + r!o/2) A/T] 

X \A{^ - 17o/2)A(f7 + f7o/2) - B{Vl - na/2)B{n + f7o/2)] . (3) 

The quantity fio — ^^o/A is the phonon energy ujq divided by the effective gap A = A(iJ) . For T ^ A, the function 
// in the integrand of Eq.(^ is approximately one in the range from J7 = to rJo/2 and zero above r2o/2. In Fig. 3 
we show some examples of the function g{ilo) for parameter values 9 = tt/2, 6 — 0.1, and A =0.1, 0.2, 0.3, and 0.6 
corresponding to H/Hc2 =0.93, 0.78, 0.61, and 0.28. For increasing A, or decreasing field, the function g tends to 
the BCS step function which is zero in the range from fio = to 2, and tt/2 for fio > 2 (see Fig. 1 of Ref. [^). The 
physical meaning of g at low temperatures for the BCS superconductor is that the minimum energy of a phonon for 
creating a pair of quasiparticles is ujq — 2A. Fig. 3 shows that, for increasing field, the relaxation rate is more and 
more smeared out and the ratio Tn/r, tends to one. A similar effect occurs for increasing impurity scattering rate 5. 

The expression for the phonon thermal conductivity limited by electron scattering and several other scattering 
processes is given by 

Kph^At^ J dxj^ — ^[l + at'^x'^ + l3t^x^ + 5tx + "ftxg{xT/A)Y^. (4) 

Here x — loq/T, t = T/Tc, and the coefficients A, a, /?, S, and 7 refer to scattering by sample boundaries, point 
defects, sheet-like faults, dislocations, and quasiparticles, respectively. For t ^ 1 and T <C A the argument of g{fla) 
in Eq.(^ is approximately equal to zero. Taking the limit fio ^ in Eq.(^ and noting that the coherence term B 
vanishes at Qq = 0, we obtain g(f^o = 0) = [A(f^o = 0, A, 9)]'^ . For constant 9 = 7r/2 this expression yields the field 
dependence of Tn/Tg for the case H || c and a cylindrical Fermi surface. For a three dimensional Fermi surface, which 
is presumably more appropriate for the quasiparticles associated with the smaller gap in MgB2, we have to multiply 
[A{0, 6')]2 by sin6' and integrate over 9 from to 7r/2. In Fig. 4 we have plotted our results for g aX 9 = 7r/2 and for 
the angular average of g vs H/Hc2 for impurity scattering rates (5=0.1 and 0.5. 

The measured Kph at H = has been analyzed in terms of an expression which is equivalent to our Eq.(^ , 
apart from the term proportional to /3. By inserting two different contributions for the electron scattering term 7 
corresponding to different gaps Ai and A2, the fitting procedure indicated that the dominant part of the scattering 
of phonons by electrons is provided by that part of the electronic excitation spectrum experiencing the smaller gap, 
A2. According to the two-band theory of Ref. [Q, this smaller gap is suppressed by a relatively small effective 
upper critical field H^2^ which leads to the measured fast drop of Kph{H) as H increases from to iJ^2^ §. We 
conclude that the measured field dependence of Kph (H) below Hj,^' at low temperatures can be explained by inserting 
in Eq.(|^) the values of the various constants obtained for i? = and, for r„/Ts = g{H), the angular average of 
5(f2o = 0) = [A(f2o = 0, A, 6')]2 . The behavior of g{H), which determines the form of the fast reduction of Kph{H), 
depends on the impurity scattering rate 5 (see Fig. 4). We note that the model of phonon scattering by vortices 
consisting of normal state cores yields a linear dependence of g on H / Hc2 ■ 

Our results have been derived for an isotropic s-wave gap. Recently an anisotropic s-wave gap A(p)/A = (1 -|- 
acos^ 0)~^/2 ];]^a,s been proposed which can explain the anisotropy of iJc2 in MgB2 [|l^. For this gap with reasonable 
parameter values of a = 10 and 20 we find that the effect of this gap anisotropy on the field dependence of C and 
Kes/'^en is similar to that of impurity scattering. For instance, the curves of C and Kes/^en versus H/Hc2 for 5 = 0.1 
and a = 20 appear to yield better fits to the experimental points in Figs. 1 and 2 than the dashed curves for a = 
and 5 = 0.5 . 

The other important effect on the upward curvature of Hes/i^en arises from the decrease of the Andreev 
scatering rate 7^ with increasing field. This decrease is given approximately by the expression 74 = 
-v/7r(A2/ sin0)exp[— (2riA/ sin0)2] where A^ ^ (-ffc2 ~ H)/H . It is interesting that the exponent in this expres- 
sion, in particular the term (sinf?)"^ ^ ig similar to the exponent in the scattering rate l/r„ which has been derived 
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previously [Q. In that model a quasiparticle is converted into a hole due to Andreev reflection by vortex screening 
currents. This process corresponds just to the imaginary part of the self energy calculated in Ref. | [lO|] . 

Our theory is strictly valid only for fields near the upper critical field Hc2 ■ However, recent numerical solutions of 
the quasiclassical Eilenberger equations for the Abrikosov vortex lattice have shown that the analytical expressions 
for the density of states ||l^,|l2| (here the function A{Vl)) yield at low energies impressively good results over the whole 
range of fields down to Hd ■ These results for low energies are important for the calculation of thermodynamic 
quantities. We expect that this is also true for the function 5(51), the spectral function of the anomalous Green's 
function. 

In conclusion, we have employed microscopic theories for the electronic specific heat C{H) and the thermal con- 
ductivity KeiH) and developed a microscopic theory for the phononic thermal conductivity Kph{H) to explain the 
measured field dependence of these quantities in the mixed state of MgB2 . We find that impurity and Andreev scat- 
tering, as well as Fermi surface and gap anisotropy, give rise to large effects. For sufficiently large impurity scattering 
rate, or gap anisotropy, one can approximately fit the data for C{H) and Ke(H) for MgB2 over a broad field range 
below Hc2 ■ Our new theory for Kph{H) is capable of describing the observed rapid reduction of the total k of MgB2 
for small fields if this is actually caused by suppression of a second smaller energy gap. 

We thank T. Dahm for helpful discussions. 
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FIG. 1. A{Q = 0, H) = js{H)/jn vs H/Hc2 for impurity scattering rates <5 = F/Ao =0.1, 0.2, and 0.5 (from bottom to top). 
Solid curves for angle = n/2 corresponding to H || c for a cylindrical Fermi surface {9 is the angle between the quasiparticle 
momentum p and the field H). Dashed curves for the average of A over corresponding to a spherical Fermi surface. The dots 
are reduced data from Ref. [8] for MgB2 at T = 3K . 
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Fig. 2 




FIG. 2. Electronic thermal conductivity ratio Kes/i^en at low temperatures vs H/Hc2 for impurity scattering rates S=0.1, 
0.2, and 0.5 (from bottom to top). Notation as in Fig.l: Solid curves for angle 6 = 7r/2. Dcished curves for the average over 9. 
The dots (triangles) are reduced data from Ref. [6] for MgB2 at T = 7.MK (T = 5.13K) . 
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Fig. 3 




FIG. 3. Ratio of phonon relaxation times Tn/rs — gi^lo) due to electron scattering at low temperatures vs Qo = cuo/A where 
uJo is the phonon energy and A = A{H) the effective energy gap. From top to bottom, A=0.1, 0.2, 0.3, and 0.6, corresponding 
to H/Hc2= 0.93, 0.78, 0.61, and 0.28. Also, 6 = 0.1 and 6 = ■k/2 (notation of Fig.l). The dashed curve refers to A = 0.6 and 
impurity scattering rate 6 = 0.5 . 
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Fig. 4 




FIG. 4. Ratio Tn/Ts = g(0,Q — 0,H) vs H/Hc2 for 5 = 0.1 (lower curves) and 5 = 0.5 (upper curves) with the notation of 
Fig.3. This quantity determines the reduction of the phonon thermal conductivity at low temperatures due to the increase of 
scattering by quasiparticles. The solid curves are the average of g over 6 corresponding to a spherical Fermi surface and the 
dashed curves are g for constant 9 = -n /2. 
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